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ABSTRACT 


Open-coast storm-surge computations are of value in planning 
and constructing engineering works in coastal regions. A quasi-two- 
dimensional numerical model for such computations is discussed from the 
standpoint of the model's underlying assumptions, range of validity, 
calibration, and application. When using simple numerical schemes 
of this sort, it is possible to make computations manually, although 
electronic digital calculations are generally preferred. 


Elementary aspects of hurricanes and the physical factors govern- 
ing storm generation processes are discussed. To display the principle 
characteristics of the model from a physical as well as a mathematical 
point of view, the basic hydrodynamic equations are given, together with 
the assumptions generally taken in their development. The equations 
consistent with the model described here are reduced forms of the basic 
equations in which several terms have been neglected. These omissions 
are indicated, and their effects on the resulting numerical scheme are 
discussed. 


The use of design hurricanes for engineering studies is treated. 
Effects of astronomical tide, initial water level, and atmospheric- 
pressure setup are considered. 


An open-coast storm-surge problem is solved for the Chesapeake Bay 
Entrance near Norfolk, Virginia. Calculations are made both by computer 
and manually. The computer program used is listed. 


FOREWORD 


This report is published to provide a simplified method for compu- 
tations of open-coast storm-surge to be used by engineers as a basis for 
construction of protective coastal works. Examples of electronic digital 
and manual calculations were used to arrive at computations based on a 
quasi-two-dimensional numerical model. 


This report was prepared by B. R. Bodine, an engineer with the 
Coastal Engineering Research Center, under the general supervision 
of George M. Watts, Chief, Engineering Development Division and 
R. A. Jachowski, Chief, Design Branch. 


At the time of publication, Lieutenant Colonel Edward M. Willis 
was Director of CERC; Thorndike Saville, Jr. was Technical Director. 


NOTE: Comments on this publication are invited. Discussion will be 
published in the next issue of the CERC Bulletin. 


This report is published under authority of Public Law 166, 79th 
Congress, approved July 31, 1945, as supplemented by Public Law 172, 
88th Congress, approved November 7, 1963. 
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STORM SURGE ON THE OPEN COAST: FUNDAMENTALS 
AND SIMPLIFIED PREDICTION 


by 


B. R. BODINE 


Section I. INTRODUCTION 


€oastal engineers concerned with design problems frequently are 
required to estimate the storm surge - the rise of water levels on the 
open coast caused by high winds acting over the Continental Shelf. On 
the East and Gulf coasts of the United States the most significant rise 
at the shore is generally associated with the fully developed hurricanes. 
Here we shall be concerned principally with those hurricanes, although 
any lesser storm is equally applicable, and surge can readily be estimated 
by the same methods. The total rise at the shore is dependent on the 
interactions of the meteorological storm with the sea and the state of 
the sea during the passage of the storm. Estimating the response of the 
sea from forces induced by the moving hurricane is complex; practical 
results are only obtained by accepting approximations. Several open- 
coast, storm-surge prediction schemes have evolved, particularly in the 
past three decades, based on various approaches, such as empirical re- 
lations, method of characteristics, statistical and numerical schemes. 
Quantitative agreement between any of the methods and observations is not 
yet completely satisfactory from the standpoint of coastal engineering 
design. In the past decade, numerical schemes have generally been used 
exclusively for predicting storm surge. 


_ High-speed data processing has brought more sophisticated numerical 
methods to the forefront in many scientific fields. Such methods yield a 
solution for the more complete and more complex non-linear partial differ- 
ential equations appropriate to the particular problem. Such techniques 
have accelerated the study of long-wave motion, and in particular tidal 
computations. Various schemes, using the method of finite differences, 
have evolved for resolving the governing hydrodynamic equations in two- 
dimensions and some in quasi-three dimensions. The open-coast surge 
problem has been treated among others by Miyazaki (1963), Leendertse 
(1967), and Jelesnianski (1966, 1967, and 1970). On a smaller scale, 
Platzman (1958) has used such a method to compute the-surge on Lake 
Michigan resulting from a moving pressure front, and Reid and Bodine 
(1968) have treated the hurricane-surge problem in Galveston Bay. Such 
methods describe the storm generation processes in a much more satisfac- 
tory manner than those more frequently used by coastal engineers. Such 
numerical schemes are mentioned here to emphasize their importance, and 
it is anticipated that in the near future they will replace the less 
descriptive schemes such as the one which will be discussed here. 


The numerical method covered here is based on the theory developed by 
Freeman, Baer and Jung (1957) which was referred to as the Bathystrophic 


Storm Tide Theory. It can be described as a quasi-two-dimensional 
numerical scheme for predicting open-coast storm surge. The method is a 
steady-state integration of the wind stress from the edge of the Conti- 
nental Shelf, taking into account some of the effects of the alongshore 
flow caused by the earth's rotation. The development of this scheme was 
an improvement over the older one-dimensional schemes which neglected the 
flow alongshore. The simplicity of the Bathystrophic approximation is 
due to the fact that integration steps can be carried out in a manner 
analogous to the one-dimensional problem. The advantage of employing 
such a scheme is that all computations can be made manually without too 
much mathematical endeavor. The computational procedure can be further 
simplified by using electronic data processing. Furthermore, the elec- 
tronic computer used does not require a very large memory capacity. 


Bretschneider and Collins. (1963) used the theory developed by Freeman, 
Baer and Jung to predict the open-coast surge at Corpus Christi, Texas, 
and vicinity. Marinos and Woodward (1968),with some modifications of the 
model used by Bretschneider and Collins, computed surge hydrographs for 
various locations along the entire Texas coast. The numerical model 
described herein is somewhat different than those used by the above 
investigators. 


There are, however, disadvantages to using simplified schemes such 
as the Bathystrophic approximation since some physical processes are 
neglected in the storm-generation problem. For some actual conditions 
this simplified scheme gives a reasonable estimate of the open-coast 
surge; for other conditions, it may not be satisfactory at all. Thus, 
at the expense of taking a more simple approach to resolving the open- 
coast storm surge problem, a computational procedure is obtained which 
may not always be useful in practice. Consequently, such a procedure is 
restricted to a certain type of surge problem, depending upon the charac- 
teristics of the shelf involved and the behavior of the meteorological 
storm. To display the limitations of the Bathystrophic theory we will 
discuss the principal underlying assumptions that lead to the theory and 
the effects resulting from the various approximations. 


Procedures and methods used to compute the open-coast surge for 
actual hurricanes and hypothetical hurricanes are covered. Also covered 
are thecriteria for establishment of the design water level at a particu- 
lar site on the coast. A numerical example of computing the open-coast 
surge is included to further illustrate the method. 


Section II. ON PREDICTION OF STORM SURGE 


1. General 


When a basin of water, such as a sea, lake, bay or estuary is dis- 
turbed by some forcing mechanism, there will be a transient motion to at 
least a portion of the water contained in the system. This sort of fluid 
motion is referred to as a surge, and can produce large variations of the 
water surface, especially in the shallower regions of the basin. Among 
the known mechanisms that cause surges are submarine earthquakes which 
generate tsunamis, a large mass of earth sliding into the water, edge 
waves, and storms. Here we shall be concerned with tropical storms, and 
generally with emphasis placed on the surges associated with a hurricane. 


The name, hurricane, applies to migratory cyclones which originate 
in the tropical regions over the oceans near the equatorial zone. Such 
- storms are referred to by different names throughout the world, although 
their origin, structure, and behavior are essentially the same. A hurrt- 
cane is called a typhoon in much of the Pacific Ocean, a tropical cyclone 
in the Bay of Bengal and the Arabian Sea, and a willy-willy off the coasts 
of Australia. A tropical storm is said to become a fully developed hur- 
ricane when the maximum wind speed equals 75 miles per hour (65 knots) 
or more than 12 on the Beaufort scale of wind force. A hurricane is 
essentially a heat engine, thus implying the existence of a heat source 
and a cold source. Its development, intensification and dissipation has 
yet to be fully understood, although considerable knowledge has been 
gained in our present century. It is beyond the scope of this paper to 
present the known complex.mechanisms which govern the hurricane; however, 
we will mention some of the elementary characteristics primarily from 
the standpoint of those related to the induced surge. 


Because of the direction of earth's rotation, the hurricane has a 
cyclonic wind circulation which is counterclockwise in the Northern 
Hemisphere and clockwise in the Southern Hemisphere. The wind system 
just above the water surface is circular, having an average diameter 
of about 400 miles across to the gale winds. of 40 miles per hour and 
approximately 100-mile diameter to the hurricane force winds. Moreover, 
the wind system extends upward in a dome-like configuration to a maximum 
elevation of about 6 miles or more. 


At the center of wind rotation and just above the surface is a 
region of relative calm called the eye of the hurricane which extends 
vertically upward through the wind system dome. The eye is cicular or 
elliptical depending on the characteristics of the hurricane, has an 
average diameter of about 14 miles, but may vary from about 4 miles to 
25 miles (Dunn and Miller, 1964). Atmospheric pressure near the sur- 
face of the ocean is minimum in the eye and usually increases in all 
directions to the storm's periphery. The temperature in the eye is 
higher than its surroundings, principally at lower levels, which causes 
the air to rise vertically, thus allowing air from the adjacent wind 


system to flow into the eye. The warm, moist air flowing into the eye 

is then heated, ascends, and is dispersed in the upper atmosphere thus 
providing a vital step in the process of transforming heat to kinetic 
energy. The more air displaced in the column of the eye, the higher will 
be the wind speeds. Since the highest pressure gradient is in a region 
near the eye, the highest wind speeds occur at the place of maximum 
pressure drop. More specifically, the maximum winds are produced at a 
region near and to the right of the hurricane center, facing the direc- 
tion which the eye is moving. 


The creation and maintenance of a hurricane requires a delicate 
balance of all processes; a breakdown in any of these processes will 
cause the storm to degenerate. The principal source of energy of the 
tropical cyclone is the release of the latent heat of condensation in 
the core or center column; thus, the hurricane can only be maintained 
over the warm ocean waters. As hurricanes move overland, they lose their 
cyclonic characteristics rapidly and transform into less intense storms 
which may move over distances of several hundred miles before completely 
dissipating. A transformation of wind speeds and direction takes place 
when a hurricane moves into a coastal area due to frictional drag over- 
land. This will be seen later to have an effect on the water level 
variation at shore. 


According to Dunn and Miller (1964), the East and Gulf Coasts of 
the United States can expect a hurricane to cross the coast on the 
average of about 2 to 3 per year. Any reaches along these coastlines 
are susceptible to tropical cyclone crossings. On the West Coast, only 
California can expect a hurricane, and only infrequently. About six 
tropical storms per year develop off the west coast of Central America 
or Mexico and only about one-half of those become full-scale hurricanes. 
Only about one in 50 of these hurricanes hits California. Tropical 
storms spawned in the North Atlantic Ocean, including the Gulf of Mexico 
and the Caribbean Sea, on the average, number about eight per year, and 
only 58 percent of these reach hurricane intensity. 


Three characteristics of hurricanes are important in establishing 
their capacity for producing storm-induced wind tides and surface waves, 
and these may be identified as: pressure differential (Ap), the differ- 
ence in pressure at the storm periphery (pn) from that at po, the cen- 
tral pressure (i.e. in the eye); radius of maximum winds (R), the dis- 
tance from the eye to the region of maximum winds; and the forward speed 
(Vp), the speed at which the eye moves. The central pressure index is 
commonly used for tagging the hurricane intensity, and is abbreviated 
as CPI. All of these characteristics appear to be dependent to some 
degree upon the storm's latitude, and they usually vary throughout the 
storm's life span. 


2. Factors Influencing the Maximum Surge 
The high winds of a hurricane which moves over the Continental 


Shelf and crosses the coast not only generate high waves but drag a 
large quantity of water along in the process. For a relatively straight 
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coastline the water transported shoreward causes the water level to rise 
both to the left and right of landfall of the storm center. Transport in 
the shoreward direction occurs in the right quandrants of the storm, and 
winds in the leading left quadrant set up a current to the left alongshore 
prior to the eye advancing into shore thus giving a rise in water level 
to the left as well as the right. For the relatively straight coastline 
the offshore winds in the left quadrants will seldom provide sufficient 
offshore transport to depress the water level below the initial level 
because of the alongshore currents developed in the leading left quadrant 
of the storm. However, if a landmass juts .out into the sea such as the 
southern part of Florida, and a hurricane center moves into the southern 
tip, then the storm driven currents alongshore become discontinuous at 
the leeward side of the landmass. Moreover when the winds in the leading 
quadrants of the storm pass over the landmass, they have a seaward direc- 
tion on the leeward coast side, and these offshore winds can depress the 
water level by several feet. All levels of surge along the coast during 
the passage of a hurricane are of concern. But from the aspect of coastal 
engineering design, it is of primary interest to know the peak surge or 
the highest surge. For the relatively straight coastline, the peak surge 
is generally assumed to occur in the vicinity of landfall of the region 
of maximum winds. This has never been fully verified from actual obser- 
vations because of inadequacy of water level recorders on the open coast, 
but merely based on the physical over-water structure of the hurricane. 
Jelesnianski (1967) has numerically shown that for a model coast of uni- 
form and straight contours, the location of the peak surge relative to 
the storm system is dependent on the angle at which the storm approaches 
the coast. Irregularities in shelf and coastline configurations would 
also govern the exact location of the peak surge. 


The highest level to which the water will rise on the open coast 
during the passage of a hurricane is dependent upon several factors or 
causes. In general, these factors can be attributed to the initial state 
of the sea, the properties of the seawater over which the storm runs, the 
behavior and characteristics of the storm, the characteristics of the 
basin and forces which act external to the system. The initial sea state 
refers to the level of the sea on the shelf heralding the approach of a 
hurricane. Generally, sea level is above normal before the arrival of a 
hurricane, especially along the coasts adjacent to the Gulf of Mexico 
and lower latitudes of East Coast. Such abnormal levels have yet to be 
fully understood, but commonly are called "forerunners", and are attrib- 
uted to transport caused by winds in advance of the hurricane. Harris 
(1963) indicates that the initial rise may be related to short-period 
anomalies in the mean sea level which are not related to the hurricane. 
In any event, whatever the cause, the initial water level can be an im- 
portant factor when attempting to account for the total rise at the open 
coast. 


A hurricane moving over the Continental Shelf affects the buildup 
of water on the coast in accordance with its intensity, path, pressure 
effect, forward speed, size, and associated rainfall. The greater the 
onshore wind speeds and usually the longer the winds act over the water, 


the greater will be the rise of water at the coast. The forward speed, 
however, can be an important factor in the total rise. If there is a 
gradual pressure change, the water level will rise in regions of low 
pressure and fall in regions of high pressure, thus the pressure is one 
of the important components of the surge. The size of the storm relative 
to the size of the basin is also an important factor, since a storm too 
large for the basin could limit the hurricane's growth as well as in- 
tensity because of the interference of the outskirt winds with the land. 
An example of a very large storm relative to the basin size was Hurricane 
Carla (1961) which made landfall near Port O'Connor, Texas. When this 
storm was approximately centered in the western end of the Gulf of Mexico 
basin, its winds at the edge were sweeping the coastal terrain of the 
east Coast of Mexico as well as the Coasts of Texas and Louisiana. 


Rainfall associated with the hurricane can contribute to the rise 
of surge on the open coast because of the volume added to the system; 
however, since the sea basin is such a large reservoir, the rise caused 
by precipitation is generally neglected. For semienclosed bays and 
estuaries, rainfall can be an extremely important factor in estimating 
the surge for any location within the system. In fact, in a semienclosed 
basin, the peak surge height may be increased as much as 3 feet due to 
direct rainfall on the water surface when coupled with the surface rain- 
fall runoff from the watersheds located adjacent to the basin. Although 
the direct effect of rainfall on the open sea can be usually neglected 
in estimating the maximum storm surge along the open coast, streams dis- 
charging large quantities of water in the sea may affect the maximum 
water level locally. 


When wind-generated surface waves of a hurricane move into the near- 
shore regions and break on the sloping beach more or less parallel to the 
depth contours, they may be responsible for a significant water transport 
shoreward. As the water moves back seaward by gravity, the momentum of 
the water particles is substantially decreased, resulting in a water sur- 
face gradient extending from the breaker region to the shore. Thus, the 
kinetic energy of the particles is increased as the wave breaks and moves 
shoreward, and decreases when the particles move back seaward. Part of 
the kinetic energy is transformed to potential energy and some energy 
is lost in the form of friction and turbulence. However, the gain in 
potential energy results in an increased water level at the shore. Wave 
setup is the superelevation of the water surface over normal surge ele- 
vation due to onshore mass transport of the water by wave action alone 
(Saville, 1962). 


Model studies of wave setup have been made by Fairchild (1958) and 
Saville (1962); theoretical studies of wave setup were made by Dorrestein 
(1962), Fortak (1962) and Longuet-Higgins and Stewart (1962, 1963, 1964). 
All of these studies show that. the water level is depressed below the 
stillwater level (SWL) at the point of maximum wave amplitude or where 
the wave peaks up and breaks. Saville (1962) indicates that waves break- 
ing on a beach face, with the crests parallel to the depth contours, 
produce a wave setup at the shore of about 10 to 15 percent of the 


incident wave height for seabed profiles 1 on 15 or gentler. Waves break- 
ing at an oblique angle to the depth contours probably produce less setup 
at the shore. Wave setup is also dependent on the configuration of the 
foreshore and the configuration of the coastline. An estimate of the in- 
fluence on wave setup of such irregularities is difficult, and would vary 
for each different configuration. 


The maximum storm-induced surge at the coast is also affected by 
the fluid properties over the total depth beneath the storm system. 
Fluid stratification is present due to the vertical density gradients, 
and strong currents in the nearshore regions cause bed materials to be 
propelled into suspension thus increasing the fluid density. Since water 
is a real fluid, there is lateral friction between the fluid stream lines. 
Such complexities are generally ignored in estimates of the maximum open- 
coast surges, although accounted for in some overall gross manner. How- 
ever, the properties of the fluid are important in the storm generation 
process, and it is not always completely satisfactory to disregard the 
various modes of internal motion caused by the disturbance of real fluid. 


The height of the hurricane surge is also governed by the charac- 
teristics of the basin, such as the size, coastline and bed formation, 
and the roughness of the Continental Shelf. As mentioned previously, a 
basin too small for the storm system can cause a reduction in intensity 
while the storm is positioned far out at sea, thereby reducing the surge 
at the open coast. The water level at the shore is also dependent upon 
the hydrographic configurations of the shelf and the shape of the coast- 
line. Thus, the water level is affected by the presence of submerged 
barriers, submarine canyons and any other formations which may abet or 
impede the flow. Normally, there is an amplification of the storm- 
induced surge wave if the shelf is relatively flat or shallow. However, 
this may not always be the case, since bottom-friction stresses may 
predominate over the surface stresses. Bottom stresses are important 
principally in shallow water and are dependent on the roughness of the 
seabed in respect to the size of the bottom material as well as the 
configuration of the bed. When estimating the surge which moves over 
the low-lying terrain adjacent to the sea, the bottom stresses become 
extremely important due to vegetation and other obstructions. 


The maximum surge on the open coast is also affected by the astro- 
nomical tide. The degree of effect depends on the amplitude of the 
tide and the phasing of the tide with the storm-induced surge. When 
maximum storm surge at the coast coincides with a maximum tide, the 
highest water level will be produced at the shore. The water level 
will be reduced when a low tide occurs simultaneously with the maximum 
surge. The vertical departure of the water surface from its mean 
position for any particular location because of tides depends on geo- 
graphical location, basin configuration, and the magnitude of the tidal 
forces at a particular time. The tide range can be as great as 18.2 feet 
at Eastport, Maine, while only 1.3 feet at Key West, Florida. Generally, 


the tides along the Gulf Coast are diurnal while those on the East Coast 
are semidiurnal. For some locations, the tide may be the principal 
component of the total rise at shore during a hurricane. 


Coriolis effect also plays an important role in the fluid motion 
beneath the storm system. Coriolis acceleration exists only when a mass 
has velocity relative to the earth's surface. The effect acts to the 
right in the Northern Hemisphere and:to the left in the Southern Hemi- 
sphere. For large basins such as the ocean, and for fluid motion propa- 
gating at relatively slow speeds, the Coriolis acceleration becomes quite 
important. 


Section III. THEORETICAL ASPECTS 


Bathystrophic Storm-Tide theory can be described as a quasi-two- 
dimensional method for evaluating the change of water level along a 
single traverse line over the Continental Shelf. Unlike the older one- 
dimensional schemes, it takes into account some of the effects of the 
earth's rotation. However, Bathystrophic theory is only an approxima- 
tion to the complete storm-generation process, and its usefulness in 
predicting actual storm surge depends upon the problem being considered. 
For some situations, the Bathystrophic approximation appears to give a 
reasonable estimate of the open-coast surge, but for others the estimate 
could be in error by a factor of two or more. Usefulness of this theory 
requires some knowledge of its underlying assumptions and the hydrodynamic 
-processes neglected in its development. To display the underlying prin- 
ciples of the approximation, the initial conditions imposed and basic 
hydrodynamic relations appropriate to the problem follow. 


nike Preliminary Conditions Imposed 


Although the basic equations which govern storm surge generation will 
not be derived here, it is important to mention the conditions imposed at 
the outset of their development. A knowledge of these conditions gives a 
better understanding of the resulting equations. The conditions are: 


a. It is assumed for a disturbance of intermediate horizontal scale 
that L >> D and L << Rg, where L is the wave length; D is the fluid depth; 
and Rp is the radius of the earth. This also implies that: 

(1) The vertical components of velocity and acceleration can 
be neglected. Thus the vertical pressure gradient. is 
hydrostatic, and vertical Coriolis effects can be ignored. 

(2) The curvature of the earth can be neglected. 


b. The acceleration due to earth's rotation is a constant. 


c. The fluid is homogenous and incompressible, thus the water 
density p is a constant. 


d. The fluid is inviscid. Thus, internal forces due to viscosity, 
are neglected. 


e. The seabed is regarded as fixed and impermeable. 


f. The effects of surface waves are considered linearly super- 
imposable on the storm surge. 


Thus, it is assumed that only horizontal flow takes place, the 
traditional approach when dealing with this type of fluid motion. Such 
flows have been referred to as nearly horizontal flows (Birkhoff, 1960). 
In the context of wave motions, such wave motions are often referred to 
as waves of long period or simply long waves. Wave motions in which the 
vertical accelerations have a marked influence on the wave behavior are 
called short-period waves or gravity waves. Surface waves referred to 
in item f above are of this type. 


The water motion which accompanies the propagation of long waves is 
unsteady, and is in a continuous state of change. This change, however, 
is not abrupt, and the motion can be considered as gradually changing. 


2. Basic Notation 


Figure 1 shows the various notations used in the discussions. The 
value D is the total depth of the fluid at time t. Moreover, D=S +d 
where S is the disturbance height of the free surface, and d is the depth 
of the undisturbed fluid. The Cartesian axes x and y are situated in a 
horizontal reference plane at the undisturbed water level, with z directed 
vertically upwards. The x and y axes are chosen counterclockwise with x 
directed shoreward and perpendicular to the shoreline. 


3. The Differential Equations 


The hydrodynamic equations may be written in either of the two 
equivalent forms: (1) the mean current velocities and (2) the volume 
transport. Preference for the particular form depends on the individual 
investigator, but generally the second form is preferred when an elec- 
tronic computer carries out the computations. Here the volume-transport 
form is taken. This form is obtained by integrating the governing equa- 
tions in the vertical over the total depth. Integrations of the primitive 
equations appropriate for the storm surge problem have been shown by 
Haurwitz (1951), Weylander (1961), Fortak (1962), Platzman (1963), Reid 
(1964), and Harris (1967). These derivations show more clearly the actual 
approximations involved. Here, the equations are taken directly in 
integrated form, since it is the purpose co display only the principal 
approximations taken. 


The governing two-dimensional equations are: 
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Figure 1. Notation and Reference Frame 
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The symbols used are defined as follows: 


U, V =x and y components, respectively, of volume transport per unit 
width. 
t = time 
Weg Wn he = Momentum transport quantities 
MO yy? XY 
fy = 2 sin = Coriolis parameter 


W = angular velocity of earth = 27/24 radians per hour: 
= geographical latitude 
tees Toy = x and y components of surface wind stress. 
Thx? "by = x and y components of bottom stress 
ry = density of water. 
Ws th = x and y components of wind speed. 
& = Atmospheric pressure deficit in head of water. 
Wy W = x and y components, respectively, of current velocity. 
\ 
P = Precipitation rate (depth/time). 
g = gravity 


Equations (1) and (2) are the expressions for conservation of momen- 
tum and often referred to as the equations of motion. These equations are 


based on Newton's law, force = mass times acceleration. Equation (3) 
is the expression for conservation of mass for an incompressible fluid, 
generally referred to as the equation of continuity. 


4. The Bathystrophic Approximation 


The assumptions made by Freeman, Baer and Jung (1957)in the develop- 
ment of the Bathystrophic Storm Tide Theory implied that: the volume 
transport perpendicular to the shore could be neglected; the onshore winds 
give an instantaneous rise of the sea surface; advection of momentum 
(field acceleration) is negligible; the sea surface is uniform, and 
parallel to the depth contours; and precipitation can be neglected. 

These assumptions indicate that: 


Tt 

dU ~ OU bx 

3t fU, Be Oo a + 0, no onshore transport, 
My? My? Myy + 0, momentum values neglected, 
as oV 


ay 5 By, + 0, alongshore sea surface uniform, 


P + 0, Precipitation neglected. 


0F cha Barometric effects are not considered here, but will be 
Bee 1 Ohy accounted for from a separate source. This effect is 


discussed later. 


Based on these assumptions, Equations (1) and (2) reduce to 
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and the continuity relation, Equation (3), is disregarded in the Bathy- 
strophic approximation because of the assumptions taken. Thus the reduced 
equations are quasi-two-dimensional since computations are restricted to a 
single axis, the x-axis; however, the rate of change of transport along 
the y-axis is retained to account for the effects of earth's rotation. 


5. Bottom and Surface Stresses 
Relations are now introduced which describe approximately the stresses 
that occur when there is water motion in the neighborhood of the seabed 


and the interaction of the winds with the sea surface. Formulas based on 
experiments and theoretical considerations have been introduced which give 
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reasonably reliable values provided that the frictional forces at the 
surface or at the bed are not influenced by each other. However, in 
shallow water, the problem is more involved since the wind effects on 
the water will influence the vertical velocity distribution over the 
total depth. Consequently, the slope of the water surface due to the 
winds is also influenced by bottom stresses. Moreover, in the vincinity 
of the shore, the yelocity at the surface may have a different direction 
from that at the bottom. Friction models which take into account the 
combined effects of surface and bottom stresses have been proposed by 
Reid (1957) and Platzman (1963). 


Here, the simplified stress laws are used, and the influence of the 
vertical velocity distribution is neglected. Moreover, similar formulas 
are adopted for representing the surface stress and bottom stress and 
this form is given by: 


2 
Teseiow (6) 


where t is the shear stress at the boundary; y is a dimensionless resist- 
ance coefficient; p is the water density; and V is the fluid velocity. 
More specifically, the shear stress at the bottom is divided by water 
density, and consistent with the stress term required for Equation (5) 
ist 


45 
Y= Kv? (7) 


Here K and V replace y and ¥, respectively, in Equation (6), and K is the 
bed friction coefficient and V as defined previously is the y-component 
of the water velocity. The bottom-stress relation given here is a dif- 
ferent form than the one originally proposed by Freeman, Baer and Jung 
(1957) and subsequently used by many investigators. As implied by the 
above authors, the bed friction coefficient had dimensions of length to 
the minus one-third power. This is in accordance with Manning's (1890) 
relation which HD Ee that the friction coefficient should be inversely 
proportional to D /3 However, on the other hand, Prandtl-von Karman 
boundary-layer theory implies that the friction coefficient is dimension- 
less (Prandtl, 1935; von Karman, 1935). For simplicity in computation, 
the latter theory is chosen. 


Equation (7) written in transport form is 
he rue pe (8) 


For typical seabed conditions, it has been found that K generally lies in 
a range from 10-5 to 5 x 10-5. 


The surface stress due to the wind formula adopted in the form of 
Equation (6) is 


tT, = pk Ww = pk ww (9) 
where W is the wind velocity. For convenience is using standard observa- 
tions, the wind velocity is taken to be at 30 feet above the wave surface, 
based on 10-minute averages. The square of the wind velocity is written 
in Equation (9) as W|W| to retain the proper sign consistent with the 
coordinate system employed. In the example calculations presented, the 
dimensionless coefficient k is taken as a function of the wind speed as 
implied by the Van Dorn (1953) relation for wind stress. Thus, it is 
assumed that 


Ie es dK for W < W 
1 c 
Ww 2 (9a) 
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ke K, + K, Goa for W > We 


where the constants K, and K, are 1.1 x 107° and 2.5 x 10a’, respectively, 


and We is a critical wind speed taken as 14 knots (about 16 miles per hour). 


Finally, surface-stress equations for introduction into Equations 
(4) and (5) can be written as: 


= = k we cos 8 
(10) 
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where 6 is the angle between the x-axis and the wind vector. 
Equations (4) and (5) can now be written as follows: 
OSs wl 2 
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Thus, two simple differential equations are obtained which can readily be 
resolved by the method of numerical integration. However, simplifying the 
more complete hydrodynamic relations (i.e., Equations (1), (2) and (3) for 
the storm-surge problem, results in formulas which are restricted to 
certain classes of problems. 


6. Limitations of Reduced Equations 


The Bathystrophic equations reveal that the time-dependent surge can 
be determined only spatially along the x-axis or a single line, often 
referred to as the traverse line. Orientation of this line becomes very 
important for valid computations due to the restrictions brought about by 
the simplifications of the resulting formulas. Since all transport along 
the x-axis is neglected, the computational axis must be oriented perpen- 
dicular to the bed contours. The onshore wind stress is then assumed 
to give an instantaneous setup of water where the deep basin beyond the 
Continental Shelf acts as a relatively infinite reservoir for replenishing 
the volume being set up. The orientation of the traverse line for shelves 
that have contours which are straight and parallel to the coastline pre- 
sents no problems. Actual shelf formations do not have such a simple 
configuration; many are quite complex in shape. If the shelf configura- 
tions are not unusually complex, particularly in the shallow nearshore 
regions, it appears that a reasonable estimate might be obtained by 
orienting the traverse line which on the average would be perpendicular 
to the bed contours, and by giving special weight to those contours in 
the nearshore region. It is emphasized that a line along which numerical 
computations are made must always be taken straight and perpendicular on 
the average to the bed contours for valid resolution. 


There are also restrictions for placement of the traverse line with 
respect to the storm system. Use of the Bathystrophic approximation is 
invalid for estimating storm surge in the left quadrants of the hurricane, 
because of the water motions due to the offshore wind components coupled 
with the effects of earth's rotation. Furthermore, the formula is gener- 
ally not considered appropriate for traverses to the left of the region 
of maximum winds. A traverse line can be taken at the region of maximum 
winds and anywhere to the right of this position, although previous 
studies indicate that surge predictions are not too reliable if the 
traverse line is selected in the vicinity of the outskirts of the storm. 
Care is needed in computing the surge along a traverse line when a storm 
moves over the Shelf with a large crossing angle. Consequently, it is 
seen that one can only investigate the dispersion of the surge at and to 
the right of maximum winds when using the Bathystrophic formulas. 


Generally, in connection with storm-surge problems, predicting the 
height of the water surface is of primary interest; predicting the associ- 
ated volume transport or current velocities is of only casual interest. 
The simplified prediction technique covered here neglects all x-component 
transport, and only partially describes the y-component of transport. 
Thus, the Bathystrophic relations should be used only for predicting the 
height of the water surface. 


This prediction scheme is considered more appropriate for determining 
storm surge for slow moving hurricanes, because of the resulting quasi- 
steady state relation. More specifically, 3U/3t is neglected while 9V/dt 
is retained. For fast-moving hurricanes, inertial effects become important, 
and can contribute significantly to the total setup at the shore. 
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7. The Numerical Scheme 


Frequently, it is of interest to know how much of the setup can be 
attributed to the onshore effects and how much can be attributed to the 
alongshore effects. To determine these two separate and distinct effects 
in the computational scheme, Equation (11) can be written in the equiva- 
lent forms as follows: 
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Equations (12), (13), and (14) can be solved by computing in incre- 
ments of time and space provided that the increments taken are not too 
large. Furthermore, the functional relationships as given by the above 
equations are taken as continuous over the entire interval. In the 
numerical analog of the equations, Ax and At are taken as nonuniform 
spacing and time steps. This will allow taking coarse spacing, Ax, where 
the seabed is relatively flat, and fine spacing where the bed slope changes 
more rapidly near shore. Nonuniform time steps, At, allow more frequent 
storm-surge computations during the period when the maximum rise in water 
level is anticipated which is sometimes useful in developing the storm 
surge hydrograph. The discrete position, x, along the traverse line and 
the time level, t, are defined as: 


Aone ) Wd) g 
i=1 (15) 
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where x, is the distance from the shoreline to the most seaward position 
prescribed on the traverse line. The summation of Ax for all 7's up to 
and including IM, the maximum value of 7 at shore, is simply x ; thus at 
the shore x = 0. Actually, according to the coordinate system taken, x 
would be negative seaward of the shoreline; however, taking x positive 
does not effect the computation scheme, and furthermore such a choice is 
consistent with the computational schemes presently used. Normally, tg 
will be taken as zero and the time level is simply the summation of At 
for all n's up to, and including, some specified value of ” = MM. 
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For the present numerical scheme, the yalues of wind stress, seabed 
depths below the undisturbed level, and coriolis effects are presumed 
supplied at all discrete positions of i, while S, D, and V are eyaluated 
at i+ 1/2. On this basis the numerical analogs of Equations (12), (13), 
and (14) can be written as follows: 
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where the overbar signifies the spatial average of B at the specified time 
levels and A and B are defined as 
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The ordinal number n represents the previous time level and n+ 1 
represents the new time level. The values of D can be regarded as the 
total water depth at a position centered between x5 and Xe The total 


water depth at any specified position not only depends upon the mean water 
depth and storm-induced setup, but also on the effects of setup due to 
pressure, astronomical forces, and the initial sea level - provided there 
is a departure from the normal level. Thus, the total water depth at a 
time midway between two time levels is given by 


pnts. pibiaoi ed i tied RA + von a i SR + (S + S Ve 
it’s 2 e 2 x y- ith 
2M (Geo Boag © WC ¢ Cal” (19) 
where Ss, = initial setup 
Sa = setup due to astronomical forces 
S,. = atmospheric pressure setup. 


Specification or determination of these setups will be discussed in more 
detail later. The total water depth at the new time level is given by 
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In other words, the total wind setup for any discrete position along the 
traverse is the setup in that reach and the accumulation of setups in all 
reaches seaward of that reach. It should be noted that in Equations (19) 
and (20), an error is introduced each time D is evaluated since the term 
(Sx + Sy) is taken at the previous time level and should be the value at 


time level n+ for Equation (19) and at n+l for Equation (20). This is, 
of course, due to the fact that the correct values are unknown at those 
time levels and an approximation is taken. This error, however, is small 
provided that increments in time and space are not taken too large. 


For problems of this sort it is customary to assume for initial con- 
ditions that the system at t = to is at equilibrium state which infers 


that the current velocities are zero and the water surface is uniform. 
For the present case, this implies that V is zero and S is uniform for 
the system. Although for actual conditions the system seldom, if ever, 
would be in a complete state of equilibrium. After a suitable lapse of 
time, this is usually of little consequence, since eventually the response 
of the system reflects only the effects of the forcing functions. 


Calculations are initiated at the seaward boundary, and then stepped 
forward through all prescribed spatial positions to the boundary at the 
shore. The same process is repeated by commencing again at the seaward 
boundary. This procedure is continued for the entire temporal range. For 
any discrete position along the traverse line, V, at the new time level is 
first evaluated based on V at the previous time level as well as B and D 
as the average value in the domain of x + Ax and t + At. A knowledge of 
V at the new time level (Equation (18)) allows determination of x and y 
components of setup (Equations (16) and (17)). The total water level 
rise at the coast can then be evaluated by a summation of the components - 
those related directly to the meteorological storm and those unrelated - 
as follows: 


S. = S. + Sy + § % + Se Sy SP Sho oe 8 (23) 


where Sy = OHg is the wave setup at the shore due to breaking waves, Hg 
is the height of the breaking wave, and % is a coefficient which norm- 
ally ranges from about 0.1 to 0.2. The local setup or setdown Sj, is 
the deviation of the water surface. from the computed water level due to 
a local effect such as inland runoff or shape of the hydrography in the 
nearshore region. Figure 2 shows a profile of the various components 
of set over the Continental Shelf at a typical coastal region. 


Since numerical calculation of V is based on a repeating formula, 
it is subject to the round-off errors, that is, each computed value will 
influence the values which are yet to be determined. To ensure that the 
value does not exceed the maximum possible value, it is useful to employ 
the following derived relationship. Equation (12) can be written 
incrementally as 


AV = k W2sine At - K V2D72 At 
or 


K v2p-2 at = k W2sino At - AV. 


For small AV, K V2D-2 ¢ k W sin. 


Thus, the y-component of volume transport becomes 
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At the new time level, the above equation can be written as 


n+l 2.n+1 
n+l HO Eanes Ores 
oe ke (24) 


Wace ll 7K 


Thus, the absolute value of the flux must never exceed the term on the 
right-hand side of the equation. This relation should always be used 

as a check, and if it exceeds this value, the flux at the new time level, 
as an estimate, can be set equal to the value given by right-hand side 

of Equation (24). 


8. Computational Formulas 


Here, we shall give the formulas in a more compact and generally 
more usable form by absorbing invariant coefficients and values caused 
by inconsistency of units into a constant for each equation. Thus, in 
the United States, wind speeds associated with hurricanes are usually 
reported in miles per hour while distances taken from hydrographic maps 
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are given in nautical miles. To eliminate conversion of units from the 
basic data, we take as follows: 


AS_, AS_ > feet 
x y 
Ax > nautical mile 
2 
g > ft/sec 
D> feet 
A, B> may nee 


mi?/hr 


< 
+ 


f > 1/hr 
At > hr 


On this basis the pertinent equations become 
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When using the above equation, one should ensure that the units used are 
consistent with those specified above. 
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Section IV. DESIGN STORMS 
1. General 


From the coastal engineering viewpoint, it is necessary to define 
the design hurricane prior to discussing the actual application of the 
numerical scheme previously described. Design analysis of coastal 
projects by engineers must provide the degree of protection necessary 
for the purpose of the project. A high degree of protection is required 
when a coastal structure provides protection of lives. A high degree 
of protection is also required when extensive property damage can be 
caused by a hurricane. When the risks are small, it may be better to 
accept a lower degree of protection to obtain lower costs for the project. 
Thus, the selection of a design hurricane for a coastal protection proj- 
ect is dependent on the economics of providing the required degree of 
protection. 


To give the most probable degree of protection required for a 
particular study area, it has been the practice to select a hurricane 
with a given set of characteristics for the particular geographical 
location. Since the characteristics are specified, the storm is called 
a Hypothetical Hurricane or Hypo-Hurricane, USA-CERC (1966). Furthermore, 
for such a storm the characteristics are taken as invariant, and the 
storm is presumed to follow a prescribed path. Two design storms, which 
are dependent upon geographical location, have been established for 
practical application in the analysis of the design of coastal structures. 
These are referred to as the Standard Project Hurricane (SPH) and Proba- 
ble Maximum Hurricane (PMH). The establishment of these design storms 
was carried out jointly by the U. S. Weather Bureau and the U. S. Army 
Corps of Engineers. 


2. Standard Project Hurricane 


Graham and Nunn (1959) developed the criteria for defining the SPH 
through a detailed analysis of past documented hurricanes. The objective 
for developing such a design storm was to provide a standard or criterion 
which would serve as a basis in the design of coastal structures subject 
to hurricane surges. In this analysis, statistical quantities were 
compiled for the experienced hurricane characteristics for various de- 
lineated zones along the Atlantic and Gulf Coasts of the United States. 
The characteristics considered were: central pressure index (CPI), radius 
of maximum winds (R), forward speed (Vp), maximum gradient wind speed 
(Vgx) and maximum wind speed (Vx). The SPH is based on a CPI having 
an occurrence probability of once in 100 years. 


Furthermore, the radius of maximum winds and forward speed are to 
be taken in such a manner that the highest surge is produced at the 
coast being investigated. The CPI, as well as R and Vp, are representa- 
tive of positions for the various zones. More specifically, Graham and 
Nunn defined the SPH as "A hypo hurricane that is intended to represent 
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the most severe combination of hurricane parameters that is reasonably 
characteristic of a region excluding extremely rare combinations." 


The maximum gradient wind speed and maximum wind speed in the belt 
of maximum winds in miles per hour were determined by the following 
formulas: 
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where K is taken as 73 for the SPH, Pp and Po are the peripheral and 
central pressures in inches of mercury, respectively, R is the radius 

of maximum winds in nautical miles, f is Coriolis parameter in units per 
hour, and Vf is the forward speed in miles per hour. Vx is the maximum 
wind speed 30 feet above the water. 


A knowledge of these variables for a given location allows computa- 
tion and graphical construction of the storm's entire isovel field. 


3. Probable Maximum Hurricane 


In the past several years, many nuclear-powered electrical generating 
stations have been constructed or proposed. Many of these have been 
sited along the Atlantic and Gulf Coasts of the United States where they 
may become vulnerable to hurricane surge and surface waves. To ensure 
that these power plants would not be flooded by hurricane surges, includ- 
ing the wind-generated surface waves, with the exception of extremely 
rare events, the Atomic Energy Commission (AEC) established a criterion 
which would provide adequate safety for public health. AEC concluded 
that adequate safety would be provided if the plant site would not be 
flooded by the surge and surface waves associated with a Probable 
Maximum Hurricane. 


The Hydrometeorological Branch of the Weather Bureau in an Interim 
Report, HUR 7-97 (1968), developed the characteristics of the PMH which 
is substantially more severe than the SPH. The PMH was defined as "A 
hypothetical hurricane having that combination of characteristics which 
will make the most severe storm that can probably occur in the particular 
region involved. The hurricane should approach the point under study 
along a critical path and at optimum rate of movement." 


HUR 7-97 not only developed the characteristics for defining the 
PMH, but updated the experience hurricane data originally carried out 
by Graham and Nunn (1959). Development of the isovel fields is essenti- 
ally the same for the PMH as that given for the SPH. The difference is 
that P, is taken as the standard sea level pressure of 29.92 inches of 
mercury for the SPH, while P, is considered as a function of the latitude 
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for the PMH. Moreover, « is considered a function of the latitude for- 
the PMH while x for the SPH was taken as a constant. 


4. Design Storms other than SPH and PMH 


It may be advantageous in some cases to select a design storm other 
than the SPH or PMH based on the risks or economics of a particular 
coastal structure. In such circumstances, the information presented in 
HUR 7-97 can be used to obtain a different probability of occurrence of 
the CPI, and obtain hurricanes of various intensities. HUR 7-97 is 
recommended instead of Graham and Nunn's earlier investigation (1959) 
because of the updating of the basic data. It should be noted that the 
CPI for a PMH was not assigned to any particular probability of occur- 
rence, but was taken as only an extreme event because of the relatively 
small amount of data presently available. Generally, if a design storm 
is selected other than SPH and PMH, it would be a storm less severe than 
the SPH. For instance, a coastal structure such as a stone jetty which 
offers protection only for normal conditions, could for some cases be 
designed for a storm less severe than the SPH. This is because it may 
be more economical to provide occasional maintenance than to build a much 
more costly and stable structure initially. A design storm for such a 
condition would be selected according to the CPI frequency considered 
satisfactory for the particular project. The radius of maximum winds 
R and the forward speed VF would be taken as those values considered 
appropriate for that particular geographical location. 


Section V. ESTIMATING THE PRESSURE EFFECT, INITIAL WATER LEVEL 
' AND ASTRONOMICAL TIDE 


1. Pressure Setup 


The effect of atmospheric-pressure setup, S,,,, which enters in the 
total depth calculation in space and time is evaluated separate from the 
relations previously described. Although several methods have evolved 
for estimating the pressure effect, none appear to be completely satis- 
factory for all storm-surge problems. The method presented here appears 
to be satisfactory for some cases, but for other cases the approximation 
may not be very reliable. Prior to giving the pressure-setup relation, 

a discussion will be made of the total problem to clarify the limitations 
of the relation presented. 


The atmospheric pressure, p, acting on the water surface during a 
hurricane is a function of position in the storm field and time, or 
symbolically p = f(x,y,t). In general, the atmospheric pressure is 
minimum at the hurricane eye, where p = pg, and increases radially out- 
ward to the edge of the storm, or storm periphery where p = p,. Because 
of the pressure variation in the open sea, water will rise in regions of 
low pressure and fall in regions of high pressure. The higher pressure 
region depresses the sea level by transporting the water to the lower 
pressure region where it causes the sea level to rise. 
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To move water from one region to the other requires time, thus the 
process is time-dependent. If the pressure change is rapid, there may be 
insufficient time for very much transport to take place resulting in only 
a small change of the water level. If there is sufficient time, an equi- 
librium state can be reached where the pressure forces are in balance with 
the gravitational forces. In the deeper regions of the sea, equilibrium 
state occurs more rapidly than in the shallow regions near the coast. In 
the nearshore regions, transport is impeded due to bed friction; more time 
is required to transfer the water. For actual hurricanes, the pressure 
disturbance moves with the speed of the storm, and can vary with the time 
elapsed. Hurricanes which move at very high speeds may not allow suf- 
ficient time for the water surface to reflect maximum and minimum levels 
that are associated with the actual pressure difference. 


Consequently, the speed of the hurricane may play a role in the 
actual amount of water being set up. Harris (1963) indicates that the 
difference in water level between any two points in the. storm field is 
proportional to the difference in atmospheric pressure provided that the 
storm speed (Vp) is small compared to the shallow-water wave speed, VgD. 
When the storm speed is comparable to the shallow-water wave speed, the 
water level will be amplified by resonance (Harris ,1957). 


The difficulty in establishing a relation which gives the actual 
amount of water set up during a hurricane is because the value has never 
seen observed for past hurricanes. The total rise at shore, even with 
the best measurement devices, is only an estimate. Separating pressure 
setup from other components of the total rise is impossible. Thus, 
methods which give only a reasonable approximation of this effect must 
be used. One such approximation that has been used frequently in prac- 
tice was developed by Myers (1954). He based his analysis on the mean 
radial pressure distribution of 69 historical hurricanes. The pressure 
p was defined as 


-R/r 
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based on the mean values of the experienced hurricanes. Here the pressure 
p is the pressure at a radial distance, r, from the storm center, pp, is 
the central pressure, R is the radius of maximum winds, and p, is the 
peripheral pressure (theoretically at r = ©). Subtracting p, from both 
sides of the equation the relation can be written as 
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When pressure is in inches of mercury, the setup of water in feet due to 
the pressure difference based on the above equation is given by 


Sy) = 1-14 (Py - Po) (i-e7B/r) (33) 
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This relation provides a method of estimating the pressure setup on the 
traverse line. However, Equation (33) does not reflect the effect of 

the speed of the storm nor does it reflect the restrictions to flow. 

Thus, it is assumed that the sea is in a state of equilibrium when con- 
sidering only the effects of atmospheric-pressure variation. For rapid 
storms, the estimate may be much too high, and when resonance occurs, the 
estimate may be too low. The best estimate possible would appear to occur 
when the storm speed is slow compared to the shallow-water wave speed. 


2. Initial Water Level and Astronomical Tide 


Because water level is often above the normal level prior to the 
arrival of a hurricane, this condition must be taken into account when 
making open-coast surge estimates for hypothetical hurricanes. It has 
been observed for past hurricanes in the Gulf of Mexico that in the 
absence of predicted astronomical tides, the water level has been as 
much as 2 feet and more above the normal level at shore just before the 
approach of the hurricane. Initial water levels observed for hurricanes 
experienced along the East Coast of the United States are substantially 
less than those in the Gulf and generally decrease with higher latitude. 
At present there are no precise methods for prescribing the initial water 
- level that should be used for a particular location in conjunction with 
a hypothetical hurricane. It has been the practice to base this initial 
level on those previously observed in the region. This appears to be the 
most satisfactory approach at this time. On this basis, when an open- 
coast surge estimate is to be made for some specific location, all per- 
tinent storm and water-level data should be analyzed to arrive at a 
reasonable value. The final value selected should also reflect the 
probability of occurrence of the hurricane event being studied. Thus, 
for an extremely rare event such as the PMH, a more conservative (i.e., 
a higher initial water level) value should be adopted provided that the 
value selected is reasonable for the area. 


Normally, in application to hypothetical hurricanes, the astronomical 
tidal rise is taken in phase with maximum level of the storm surge. For 
an SPH, the spring high tide is usually taken; for the PMH, that tide 
exceeded by only 10 percent of the spring tides is taken. The specific 
tide amplitude chosen depends on the particular geographical location 
being studied. 


Including the effects of the initial water level and astronomical 
tide in the numerical integration procedure gives a more accurate esti- 
mate of the water-level rise at the shore. This is due to more correctly 
specifying the total water depth for all steps in time and space, and the 
errors introduced by directly adding these effects to the wind surge are 
generally very small. Thus, in most cases it is usually satisfactory 
to include these effects after making the storm wind surge calculations. 
Another method frequently employed is to include the assumed astronomical 
tide level and the initial water level as constants in numerical inte- 
gration procedure and correcting the final water level hydrograph for 
the time-dependent tide. Taking the tide as a constant for the entire 
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computational time reduces the work involved in ensuring correct phrasing. 
This scheme also involves errors in the computations due to incorrect 
depths, but for evaluation of the peak surge the solution would be more 
satisfactory. In the examples we will use the latter technique. 


Section VI. APPLICATION 
1. General 


Because the numerical model described here is quasi-two-dimensional, 
calculations can be carried out manually, although for most problems 
there is a considerable amount of work involved. Use of high-speed 
electronic computers for performing the storm-surge calculations not 
only reduces the work, but reduces the chances for error. Thus, when 
electronic computers are available, machine computations are generally 
preferred, particularly when several storm-surge estimates are to be made. 
In the case of either the SPH or PMH, several separate computations are 
generally needed to ensure proper selection of the storm characteristics 
for a specific location on the coast. Consequently, manual computations 
would be very laborious. However, the method for carrying storm-surge 
computations manually is shown by an example as well as by machine 
computations. 


The formulas presented have been discussed from the viewpoint of 
evaluating water elevations associated with hurricanes. However, nothing 
restricts the formulas from being used for predicting a surge from extra- 
tropical storms provided the storm moves in a shoreward direction over 
the Continental Shelf. For some coastal regions where hurricanes are not 
experienced or where strong onshore winds may produce higher surges than 
hurricanes, the numerical method presented may be useful in predicting 
setup at shore arising from high winds. Computations of extratropical 
wind surge in regions where water levels associated with hurricanes are 
dominant, may provide insight in calibrating the numerical model for a 
specific location. Thus, estimates of open-coast surge can be made for 
wind storms as well as hurricanes. 


2. Storm Surge Computer Program 


A program listing for estimating the open-coast storm surge is given 
in the Appendix. This program is written in FORTRAN IV, for the UNIVAC 
1108 and EXEC II operating system. The UNIVAC 1108 used has 64,000 words 
of core memory and additional memory from magnetic drums and tapes. 
Numerous comment cards have been provided in the program to identify the 
symbols and to give a running commentary of the procedures employed. 

This program with its moderate memory capacity requirements and program 
simplicity can be readily adapted to other machines in most cases. 


Input to the program is supplied by cards. Considerably more cards 
are required for computing storm-surge associated with actual or observed 
hurricanes than for computing surge associated with hypothetical hurri- 
canes such as the SPH or PMH. Much of the basic information needed for 
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estimating hypo-hurricanes is generated internal to the program. In 
subsequent discussions of an example problem, the basic input to the 
computer will be more clearly demonstrated. 


3. Example Problem 


For illustration purposes, the peak surge will be determined at the 
mouth of the Chesapeake Bay using characteristics given by Graham and 
Nunn (1959) for a storm with a large radius and a moderate forward speed. 
It is not implied here that this storm is a Standard Project Hurricane 
(SPH), since a study would be needed to determine if it produced the 
highest possible peak surge at this location. The parameters taken are 
as follows: 


CPI = 27.57 inches of mercury 

P, = 29.92 inches of mercury 

R = 35 nautical miles 

Vp = 22 knots (25.3 miles per hour) 
Vx = 102 miles per hour 


The isovel pattern considered appropriate for the above parameters is 
given by Graham and Nunn (1959) on their Figure 33. 


Figure 3 shows the Chesapeake Bay Entrance and the offshore hy- 
drography over the Continental Shelf. The line along which computations 
are carried - the traverse line - corresponds to a latitude of 37°00'. 
This traverse line represents a line which, on the average, is about per- 
pendicular to the bed contours. The path of the storm is taken parallel 
to and about 35 nautical miles south of the traverse line. Since com- 
putations of the surge would be invalid within the Bay, the last raw 
data position is at the mouth at a longitude of 76°00'. Figure 4 shows 
the approximate bed profile along the traverse line from the mouth of 
Chesapeake Bay seaward to the 600-foot depth contour. Experience has 
shown that a point where the depth is 300 feet is usually sufficient for 
commencing the computation from the seaward position to the most landward 
point. Setup in regions deeper than this depth are generally negligible. 
This will become evident by the example. The discrete points in space 
in which raw data is supplied and the corresponding depths below mean 
sea leval are taken as shown on page 53 in the Appendix. 


By positioning the isovel so that the maximum wind lies on the 
traverse line, and winds are about zero at the most landward position, 
the values of wind (W), radius (r), and theta (0) can be plotted as a 
function of distance from the mouth of the Bay seaward. For the example, 
this is shown on Figure 5 where time in hours indicates the position 
relative to the mouth after an elapse of time. The method of evaluating 
W, r and @ along the traverse line is shown on Figure 6. 
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Figure 3. Offshore Bed Contours and Traverse Line for Chesapeake Bay 
Entrance 
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Figure 6. Radial Distances and angles of the Wind Along Traverse Line 
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From Figure 5, the parameters shown can be read for all selected spatial 
points and at every time step used. In actual practice, Figure 5 should 
be plotted on large rectangular graph paper so that the specific values 
can be read more precisely. The one shown here is for illustration 
purposes only. 


Resolution of the problem will be first considered from the stand- 
point of computer applications, and then from the standpoint of manual 
calculations. 


The input data in connection with computer program for the example 
problem is shown in the Appendix on page 52. Every line of data repre- 
sents an input card, and each card or set of cards is numbered at the 
right for reference to the following explanation: 


CARD 1: This is a title card which is to be used for identification 
of the problem, such as location, date, etc. 


CARD 2: The first number indicated as 62 on the card is the total 
number of time steps used in computing the storm surge, or symbolically, 
this is denoted as NM = 62. With reference to Figure 5, it is seen that 
if 35 hours of computational time are taken, then based on the forward 
speed of the storm of 22 knots, the distance required on the graph 
for computational purposes is Vpt = 770 nautical miles. However, the 
parameters must be known 62 nautical miles beyond this distance due to 
the width of the shelf or a total distance of 832 nautical miles. If 31 
hours are selected as computing time, it is seen that the parameters are 
always defined, thus this time will be used in the computations. Further- 
more, if computations are carried out in time steps of one-half hours, 
then the required number of time increments is 62. 


The second value is the number of shelf reaches (LM). There are 18 
positions on the traverse line for supplying raw data (see page 53 of 
the Appendix), thus there are 17 increments along the line. 


The third value is an option of the manner in which data is to be 
supplied. If IDATA = 1, all wind-stress, pressure, and astronomical- 
tide data are supplied by card externally, otherwise the data are 
generated internally. 


The fourth and fifth values, IOMIT and IOMITD, if equal to 1, omits 
a detailed printout of output and input data, respectively; otherwise a 
detailed printout is made. 


CARD 3. The set of cards shown here contains the values of the time 
increments (DELT). Here, they are taken identical for all NM. However, 
the program is written so that any time increment can be taken for a 
given time level. 


CARD 4: The cards in this set give the distances in nautical miles 
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(DIST) from shore for all discrete points along the traverse line be- 
ginning with the most seaward point. 


CARD 5: All cards in this set give the depths (D) below the still- 
water level for all points corresponding to DIST above. 


CARD 6: This set represents the latitudes (PHI) at the discrete 
points along the traverse line. Since the traverse line corresponds to 
a latitude 37 degrees, all values are the same. 


CARD 7: The first five values on this card, R, Po, Pn, VF and Se, 
were defined previously. The next three values are the constants Cj, C2, 
and C3, which correspond to the constants given in Equations (25), (26), 
and (27), respectively. The last value of 1.1 is a wind-stress correc- 
tion factor (WKCOR). This increases Van Dorn's (1953) wind-stress factor 
by 10 percent to account for the additional stress caused by energy 
imparted to the sea due to precipitation. 


CARD 8: The first value is maximum gradient wind speed (WX) and the 
second value is the energy dissipation factor. Since calibration of the 
numerical model has not been carried out for this location, a value of 
0.0025 is taken. 


CARD 9: The cards in this set give the height of the astronomical 
tide above the mean sea level datum. Here the astronomical tide is taken 
as a constant for the entire computational period with anticipation of 
correcting the final hydrograph. However, the program is written in such 
a way that an actual tidal hydrograph can be provided in the input. 


The remaining cards are provided only if IDATA is not equal to l, 
which is the case for this example. These cards, with the exception of 
the first one, give the distances indicated on the horizontal axis on 
Figure 5 and the corresponding values of the parameters shown. The dis- 
crete values of r, W and © are chosen in such a way that the intermediate 
values can be interpolated by the program linearly. Thus, when the curves 
change rapidly, values must be prescribed frequently with distance; when 
curves change slowly, values can be prescribed less frequently for rela- 
tively straight curves. The program by an interpolation procedure recon- 
structs the curves in tabular form for each unit of distance for possible 
use in generation of the basic storm surge information. 


CARD 10: The values 21, 28 and 26 are the number of points selected 
for adequate representation of Tf, W and © curves, respectively. 


CARDS 11, 12, and 13: These cards give the distances along the 
horizontal axis in which r, W, and 6 are specified, respectively. 


CARDS 14, 15 and 16: These cards give the values of r, W, and 6 at 
the specified points along the traverse line, respectively. 
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The output of the program is shown on pages 54 and 55 in the 
Appendix. Page 53 shows a printout of the basic information of the open- 
coast storm-surge problem. For example purposes, the values indicated 
below the series of stars on page 53 are the wind stress, astronomical 
tide, and pressure setup data generated internally in the program for 
all reaches along the traverse line for two time levels. These data 
provide the basic information for the storm-surge problem, and is input 
by card when IDATA = 1. In case that IOMITD = 1, this information is 
not printed. 


Page 54 in the Appendix shows an example of a detailed printout of 
the storm-surge results for a single time level. If IOMIT = 1, this 
information will not be printed. 


Page 55 in the Appendix gives a summary of the problem, together 
with the various components contributing to the total setup at the 
mouth. Thus, it is seen that the peak surge for this storm for the 
Chesapeake Bay Entrance is 13.4 feet. The estimate of the peak surge, 
tabulated in the column of total water level, is assumed to be correct, 
but other values for the total water level would require correction for 
the tidal oscillation. Such an adjustment will not be shown here since 
it is only a matter of trigonometry and arithmetic. 


Because the surge is not being evaluated at the shore for this 
particular problem, wave setup is not a factor in the total setup. Other 
factors such as effects of the inlet, local bed configurations, may affect 
the total water level at the mouth, but since we don't know how to esti- 
mate them, we will consider these effects negligible. Thus, the peak 
design water level would be taken 13.4 feet for the given storm at the 
mouth of Chesapeake Bay. 


Bretschneider (1959) used a one-dimensional numerical scheme for 
computing the open-coast storm surge for the Chesapeake Bay Entrance. 
Although he used slightly different hurricane characteristics than the 
ones used here, it appears that for this particular case that both 
schemes would, for all practical purposes, give a wind surge on the 
same order of magnitude. 


We will now consider determining the storm surge in the Chesapeake 
Bay Entrance by manual calculation. For such problems, a bookkeeping 
system is needed, so that all the values to be calculated can be carried 
out in a systematic manner. One method is to construct tables similar to 
the one shown on Table 1 for each time level. Such tables allow stepping 
through all discrete points in space as well as time. For example pur- 
poses only, Table 1 shows the necessary computations at a time level 
centered between t = 16.5 hours and t + t = 17.0 hours. This corresponds 
to the time of peak surge found previously from the results of the com- 
puter program. The table not only shows the basic data for the present 
time level, but the necessary values evaluated at the previous time level. 
Many of the columns give a breakdown of the pertinent equations presented 
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in the text so that the final equations can be resolved systematically 
in steps. 


Finally, we see that the peak surge estimated manually is the same 
as that estimated by high-speed digital computations. However, there is 
usually a difference in the two values because of truncation and round-off 
errors. 
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TABLE 1 


MANUAL COMPUTATIONS AT A TIME LEVEL CENTERED BETWEEN 
t = 16.5 hours and t + At = 17.0 hours 
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5.0 60.0 0.60182 0.5 2.5 6.730 1.71 0.0384 0, 0257 32.0 42 101.5 1,094 0,335 «1.78 7656 6894 = -2.87x10°° 0, 0242 
10.0 56.0 32.0 46 101.0 2 
5.0 50.0 0.60182 0.5 2.5 7.495 1.65 0.0332 0.0263 32.5 so 100.5 1.077 0.341 «1.77 6492 7737 2.87x10°© 0.0205 
5.0 440 33.0 53 100.0 
1s 42.0 0.60182" 0.5 2.5 + 7.714 1.60°" 0.0290 0.0263 * 33.2 52 99.8° 1.053 0.349 1.74 6132 7849 2.86x10"© 0.0193 
3.5 40.0 33.5 56 99.5 
1s 38.0 0.60182 0.s 2.5: 7.929 1.58 0.0269 0.0262 33.7 57 99.2 1.037 0.355 1.73 $360 8253 -2.86x10"° 0, 0169 
2.0 36.0 34.0 58 99.0 
1.0 34.5 0,60182 0.5 2.5 8.065 1.56 0.0243 0.0247 34.4 58 96.1 1.017 0.362 1.71 4894 7832 2.84x10"° 0.0153 
1.0 33.0 34.8 59 93.1° 
0.5 31.5 0.60182 0.5 2.5 8.189 1.55 0.0123 0.0219 34.9 60 90.2 1.003 0.367 1.70 4068 “70a = 2.79x10° 00125 
0.0 30.0 35.0 60 87.20 
NOTE: 


Time (t) = 16.5 hours; Time (t+at) = 17.0 hours. 
at = 0.5 hours. 

Storm Characteristics and other Basic Information: 
P, = 27.87 inches of mercury @ = 37.0° + sing = .60182 (constant in this case) 
P,, = 29.92 inches of mercury k = 0.0025 
Vp = 22.0 knots Wind Stress Correction factor = 1.10 


R, = 35.0 nautical miles 4p = (Pp-Po) = 2-35 inches of mercury 
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Columns 6 through 9 are values evaluated at previous time level. 

Columns 10, 11 and 12 are values taken from curves given on Figure S. 

* Wind reduced by a factor G£ 0.945 due to interference with land mass. 
** Wind reduced by a factor of 0.890 due to interference with land mass. 


This is a test. 


Table 1 


(Continued) 


26 27 29 an 
20 21 22 23 24 25 aac ze 
neh nel s nel +1 
nel D D 1+ Col.25 Col.24ev"  * * re ve 
er nel Eq. (19) Eq. (20) yes 2 2 (Value 
+s Col.26 (5280) nel 
1.10: k 2) Jcois. 2, Scots. 2, kat Eq. (27) (o"**) used) 
+ Col. 17 z 4,5,6,21 4,5,6,15 7 2 
(Mi/tr)? (ft) (ét.) (ft.)  (Mi/Rr) “Wl i? /Hr) * (iC /Hr) 
0.0151 1.31 480.12- 480.04 -0.0078 120.94 1.0120 0.0710 69,696 3479.62 0.2234 0.0710 
0.0156 1.35 299.76 299.69 -0.0077 310.26 1.0290 0.0651 69,696 1401.10 0.1418 0.0651 
-0.0156 1.42 212.97 212.90 -0.0075 614.65 1.0299 0.0305 69,696 707.09 0.1007 0.0308 
-0.0148 1.49 171.83 171.76 -0.0060 944.21 11,0712 0.0507 69,696, 436.62 0.0791 0.0507 
-0,0125 1.55 148.35 148.28 -0.0031 1266.76 1.0089 0.0025 69,696 274.84 =~ -0.0628 0.0025 
-0.0088 1.63 129.64 129.57 -0.0001 1658.78 1.0939 0.0413 69,696 147.74 0.0460 0.0413 
0.0006 1.69 116.50 116.44 0.0035 2054.07 1.1363 0.0498 69,696 8.13 0.0108 0.0108 
0.0080 1.71 101.88 101.88 0.0067 2685.90 1.1763 0.0503 69,696 83.04 0.0345 0.0345 
0.0139 1.77 87.32 87.32 0.0090 3656.29 1.2130 0.0458 69,696 105.98 0.0390 0.0390 
~ 0.0184 1.77 78.19 - 78.20 0.0168 4560.01 1.2405 0, 0476 69,496 112.520.0402 0.0402 
0.0218 1.75 71.48 71.51 0.0119 $456.31 1.2619 0.0399 69,696 111.48 0.0400 0.0399 
0.0244 1.71 62.21 62.27 0.0127 7203.56 1.2989 0.0353 69,696 94,61 0.0368 0.0353 
0.0260 1.67 54.38 54.45 0.0131 9427.35 1.3417 0.0314 69,696 77.08 0.0333 0.0314 
0.0260 1.65 50.58 50.66 0.0131 10897.08 1.4364 0.0278 69,696 66.73 0.0309 0.0279 
0.0245 1.63 47.20 47.28 0.0123 12513.65 1.3801 0.0265 69,696 54.77 0.0280 0.0265 
0.0216 1.62 44.31 44.39 0.0109 © 14199.22 1.2183 0.0190 69,696 42.56 0.0247 0.0190 
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32 33 34 35 36 37 
poy pment ea Total 
ax 34nd + Col.32 = Jas ox A S = JS. Setup 
Tel a yned (21) o”* by GE SF 
Eq. (25) Eq. (25) 4 ) gq. (23) 
(£t.) (£t.) (£t.) (£t.) (£t.) 
0.442 0.038 0.038 0.029 0.029 4.30 
1.768 0.139 0.177 0.131 0.160 4.62 
2.489 0,091 0.268 0.222 0.382 5.00 
3.086 0.188 0.456 0.304 0.686 5.56 
3.574 0.011 0.467 0.386 1.072 6.02 
4.090 0.203 0.670 0.481 1.553 6.78 
4.552 0.059 0.729 0.567 2.120 7.48 
5.202 0.216 0.945 0.638 2.758 8.41 
6.070 0.285 1.230 0.695 3.452 9.45 
6.777 0.328 1.558 0.709 4.162 10.50 
7.412 0.356 1,914 0.687 4.849 11.54 
8.511 0.362 2.276 0.668 5.517 12.56 
2.920 0.110 2.386 0.216 5.733 12.86 
3.139 0.105 2.491 0.203 5.936 13.16 
2.242 0.072 2.563 0.131 6.067 13.34 
1.196 0.027 2.590 0.057 6.124 13.41 
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APPENDIX 


Storm-Surge Program Listing . 
Cardminput Data cy. sin ee mete 
Basic Storm-Surge Information . 


Optional Detailed Output 


Summary of Open-Coast Storm-Surge Problem . 


45 


Page 

46-51 
52 
53 
54 


55 


AAANDANDADAAANAADADAANADAADAANADAAANAADAADAANAAAANAAGDLGAAAAAH 


NOTATION 


x< 
R 
WW 


L nl 
cx) 
Woot We a Ot at at ta 


RS 
DELT 
Vv 

PHI 

PI 
VF 
wKCOR 
1,¢27¢3 
wX 
wK = V 

OFF = 
Or LSXN 
OrLSY 
Sx 

SY 

SPP 
SPN = 


U 
S 
WGP PP OP Pe th Pt ht eat ta 


I 
I 
1 
T 
p 
Pp 


Wo bi atte 


QUAST=TWO=DIMENSIONAL OPEN COAST STORM SURGE PROGRAM 


USED: 

REFERENCE HORIZONTAL DISTANCE FOR RADII IN NeoMe 

REFERENCE HORIZONTAL DISTANCE FOR WIND IN NeMe 

REFERENCE HORIZONTAL DISTANCE FOR THETA IN NeMe 

TOTAL NO. OF RADII VALUES READ FROM CURVE. 

TOTAL NOe OF WIND VALUES READ FROM CURVE. 

TOTAL NOe OF THETA VALUES READ FROM CURVE. 

RADII READ FROM CURVE IN NeMe 

WIND SPEED READ FROM CURVE IN MPHe 

THETA READ FROM CURVE IN DEGREES. 

ANGLE MEASURED FROM THE POSTIVE X “AXIS COUNTER CLOCKWISE 
THE WIND VECTOR? (THETA). 

DISTANCE FROM COAST IN NeMe 

TOTAL WATER DEPTH IN FEET. 

SETUP DUE TO ASTRONOMICAL FORCES IN FEET. 

INITLAL SFTUP IN FEETe 

SETUP DUE TO PRESSURE IN FEET. 

DEPTH BELOW MEAN WATER LEVEL IN FEET. 

WIND SPEED IN MPHe 

W *« w * COS(THETA) IN MPH SQ. 

Wo o« WwW * SIN(THETA) IN MPH SQ. 

hMAXe NOe OF TIME INCREMENTS. 

wAXe NOe OF SHELF REACHES. 


~CENTRAL PRESSURE INCHES OF HG. 


PERIPHERAL PRESSUE IN INCHES OF HGe 

RADIUS OF MAXe WINDS IN NoMe 

RAPII MEASURES FROM STORM CENTER TO POINT ON TRAVERSE TN 
TimE INCREMENTS IN HOURS. 

Y-COMPONENT OF VOLUME TRANSPORT. 

GEOGRAPHICAL LATITUDE IN DEGREES. 

3e14159 

FORWARL SPEED OF STORM IN KNOTS. 

CORRECTION FACTOR FOR WIND STRESS COEFFICIENT. 
DIMENSION CORRECTION FACTORS 

MAXIMUM GRADIENT WLND SPEFD (MPH) e 

AN DORNS WIND STRESS FACTOR. 

BOTTOM FRICTION FACTOR. 

NCREMENTAL X-COMP SETUP IN FEETe 

NCREMENTAL Y-COMP SETUP IN FEETe 

OTAL X=-COmP SETUP AT SHGRE IN FEET. 

OTAL Y-COmMP SETUP AT SHURE IN FEET. 

RESSURE SETUP AT PREVIOUS TIME LEVEL. 

RESSURE SETUP AT NEW TIME LEVEL. 


TO 


NMe 


DIMENSLON SP (DU) eww (50) eWwY (5G) eWK (50) 00050) eDIST(50) eXR(100) 0 Xwl 
1100) oXT(100) oRS(1500) WC L50U) eV (10U) eVW(100) oVT(100) 0A(20) »SA(G9) 
200FLT (200) oPHI (50) 0161500) UAVG(S0) pDELX (50) 65X(100) »SY(1N0) »>SWw(d 
300) eSPRES(10U) 9STOT(10G) oe STNPHI (50) 6SxXP (100) oSYP(10U) »SPP(10N) oP 


4100) 0& 
READ STOR 
READ ( 

1 FORMAT 
READ ( 

40 FORMAT 
READ ( 
REAU ( 
READ ( 
REAL ( 

45 FORMAT 


P(100) 

M ITOENTIFICATION TITLEe 

501) (All) »I=192u) 
(20A4) 

5°10) NMelLMe TOATA® TOMIT» TOMITO 
(515) 

5015) (UEIT(i)» LT=19NM) 

50i5) (DIST(L) eo L=LeLM) 

5015) (UCL) eL=1°LM) 

5015) (PHI(L) eL=le LM) 
(8F10-4) 
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LMM = LM = 1 
PI = 3.14159 
READ (5°20) RePOsPNeVFeSteCleC2.C3ewKCOR 
READ (5°2U) WXeRFF 
20 FORMAT (10F7.2) 
READ (5015) (SA(N) »N=10NM) 
Cc IF IDATA = 1 eoWIND STRESS» ASTRONOMICAL SETUP»PRESSURE SETUP AND WINN 
Cc STRESS COEFFICIENT DATA SUPPLIED BY AN EXTERNAL SOURCE» OTHERWISE 
Cc DATA GENERATED INTERNALLY FROM WINDeRADII* AND THETA CURVES. 
IF (IDATA -EQe 1) GO FO 220 
READ (5°10) IReIwoIT 
READ (5°30) (XR( I) eT=le IR) 
READ (5030) (XW(I) eI=lo IW) 
READ (5030) (XT( 1) eI=leIT) 
READ (5°30) (VR(I)»I=1e IR) 
READ (5°30) (VW(I)»I=1eIw) 
READ (5030) (VT(1) eI=l1eIT) 
30 FORMAT (1068.2) 
XRE = XR(IR) + 120 
IND = XRE 
DO 35 I=i1e1T 
35 VTI(I) = (VT(1) * PI)/180-0 


Ik = 1 
C IN THE FOLLOWING THREE LOOPS ALL ORDINATE VALUES ON THE CURVES ARE 
C INTERPOLATED FOR EACH UNIT OF DISTANCE> 


DO 60 I=10eIR 
RS(IK) = vR(1) 
{I = XR(1) 
IKN = XR(I+1) +41 
III = IKN = 1 
IF (II e£Q@. IIL) GO TO 50 
Ik = Ik + 1 
RATIO = (VR(It1L) = VR(T))/(XR(141) = XRCI)) 
uO RSCIK) = RS(IK=-1) + RATIO 
Ik = Ik + 1 
IF (IK eGée IKN) GO TO 60 
GO TO 4U 
50 IK = Ik +1 
60 CONTINUE 
IK = 1 
bDO 90 I=1eIw 
w(IK) = vw(1) 
II = Xw(1) 
IKN = XW(I4+1) + 1 
III = IKN = 1 
IF (II eEQ@. III) GO TO 80 
1K = IK + 1 
RATIO = (vw(It1) = vw(l))/(Xwlr41) = xwlt)) 
70 w(K) = wlIK=-1) + RATIO 
IK = IK + 1 
IF (IK eGE. IKN) GO TO 90 
GO TO 70 
g0 IK = IK + 1 
90 CONTINUE 
Ik = 1 
KiI = 0 
DO 120 I=isIT 
IF (KII .ta@e 1) GO TO 93 
TCIK) = vT(1) 
IF (VT(I+1) »-GT. 020) GO TO 95 
Pitse Sos 
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an 


io) 


io) 


ioe) 


93 


S5 


100 


110 
120 


IN 


GO TO 9 
VT(I) = 6.2832 
TOLER) =v Gy) 
Kees 

rey Se SOG ed) 

IKN = XT(I4+1) + 1 

IIt = IKN = 1 

IF (II ef£Q. III) GO TO 110 

Ik = Ik + 1 

RATIO = (vTCItl) = vT(1I))/OXT(141) = xTCI)?} 

T(IK) = T(IK=1) + RaTIO 

MEGS US < R 

if (IK eGee IKN) GO TO 120 

GO TO 1040 

ik = Ik + 1 

CONTINUE 

HOUR = 0.0 

THE FOLLOWING Loop THE PERTINENT INTERPOLATED VALUES ARE WSEN To 


GENERATE THE BASIC INPUT DATA FOR THE STORM SURGE PROBLEM. 


DO 210 N=leNM 


THE STORMS POSITION RELATIVE TO THE SHORELINE IS DETERMINED. 


DISTC = HOUR * VF 

DO 170 L=ieLM 

DI = DIST(L) + DiSTC + 160 
Lb = DI 


SETUP DUE TO THE ATMOSPHERIC PRESSURE OIFFERENTIAL IS EVALUATFD. 


SP(L) = 10614 * (PN = PO) * (1060 = EXP(=R/RS(LD))) 
IF (DIST(L) =-1-0) 130°1400145 


THE wIND SPEED NEAR SHORE DUE TO THE STORMS INTERFERENCE WITH THE 


136 
140 


145 


150 


LAND MASSe 

WIND = W(LOD) * 0.89 

GO TO 150 

WIND = W(LD) * 02945 

GO TO 150 

wIND = W(LD) 

wSQ = WIND * WIND 

WWX(L) = wS@ * COS(T(LD)) 
wwY¥(L) = wS@ * SIN(T(LD)) 

iF (WILD) .~LEe 1660) GO TO 160 


VAN DORNS wIND STRESS COEFFICIENT IS DFTERMINED. 


160 
170 


200 


210 
220 


wK(L) = 6-0000011 + 060000925 * (1606 = 16.0/W(LD) ) *¥*2 

GO TO i70 

WK(L) = ¥e9000011 

CONTINUE 

WRITE (39) DELT(N) 

DO 200 LL=1»LMM 

WRITE (39) WwX(LL41) oWwWX(LL) oWWY(LL+1) oWWY (LL) 0 SAC(N) oe SP(LI+1) e SPC 


1L) »WK (LL) 


CONTINUE 

HOUR = HOUR + DELT Cin) 
CONTINUE 

CONTINUE 

WRITE (60230) 


THE BASIC INFORMATION OF THE STORM PROBLEM IS WRITTEN OUT. 


230 


FORMAT (1h1030X,52HQUASI-TWO=DIMENSIONAL OPEN COAST STORM SURGE PR 


1OGRAM///) 


WRITE (601) (ACI) 0 1=1920) 

WRITE (69250) 

FORMAT (54xX010HINPUT DATA0/*?54X0910H===== m=-=9///) 

WRITE (60260) POsPNOROWX0 VE OSE RFF oe WKCOR*C19C20C3 

FORMAT (11X019HCENTRAL PRESSURE = 0F50e298H INe HGesr22Xe22HPERIPHFR 
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LAL PRESSURE = 9FSe20GH INe HGe//05Xe25HRADIUS TO MAXIMUM WIND = oF 
24e1095H NeMer24Xe e4HMAXIMUM GRADIENT WIND = oF Sele 7H MI/HReo//010X0 
S20HTRANSLATION SPEED = 0F4eleGH KNOTS*#32Xe1SHINITIAL RISE = »F4.20 
45H FEET®//e5Xe?25HBOTTOM FRICTION FACTOR = eF6.4019Xe32HWINN STRESS 
5 CORRECTION FACTOR = Fue2e//017X»5HCL = ePF6e1e21XeSHC2 = oF6.1021X 
6eSHCS = oF Tel e///) 
WRITE (60270) 
270 FORMAT (40Xe41LHSFABED PROFILE OVER THE CONTINENTAL SHELFo//e46X013 
LHDISTANCE FROM? 8X er SHDEPTH®/*(47Xe1LAHSHORE (NeMe) »8Xe6H(FEET) ) 
DO 280 L=irelLu 
WRITE (60290) DIST(L) eU(L) 
260 CONTINUE 
250 FORMAT (49xX9F5e209KeF Sel) 
IF (IDATA .NEe 1) GO TO 340 
DO 33U N=10NM 
STORM DATA IS READ FROM CARUSe THIS OPTION IS USED FOR SUPPLYING RAw 
DATA ASSOCIATED WITH ACTUAL STORMS. 
READ (59300) DELT(N) 
300 FORMAT (6xXF6é.2) 
WRITE (39) DELT(N) 
bo 320 LL=1°LMM 
READ (50310) WWX(LL+1) eWWX (LL) owWWY(LL4+1) owWWY(LL) ®@SAC(N) o SP(LL41) 0 SP 
L(LL) oe WK CLL) 
WRITE (39) WwX(LL4+1) »WWX(LL) oe WWwY (LL4+1) eWWY(LL) » SA(N) op SP(LL41) oSP(L 
1L) epwK (LL) 
310 FORMAT (4F10.006X03F602°?5X0F1009) 
320 CONTINUE 
330 CONTINUE 
340 CONTINUE 
REWIND 39 
TII = Oecd 
IF (IOMITD -EQGe 1) GO TO 405 
IF IOMITD NOT EQUAL 1e INPUT COMPUTATIONAL DATA PRINTED OUT. 
bO 400 N=leN- 
WRITE (60350) 
350 FORMAT (//7/01X%011 91% ARO i IO OI OR fk ROK IO ga ge ik a i ak ofc 2 
De RR ORK CIO IG AI I i ROK iO io i A OR Rk a i I go koi ae ie a Re ake age 
2 KK KKK ) 
READ (39) DELT(N) 
Tie Ss Wate 
TII =TI + DELT(N) 
WRITE (60360) TIeTII 
360 FORMAT (1X0/i1X11HTIME (T) = oF6e.206Xe16HTIME (T+DELT) = oF Ge26//) 
WRITE (69370) 
370 FORMAT (10X09 8HWWX(I+1) e 7XeGHWWX(1) 0 7Xe BHWWY (IT+1) oBXoGHWWY (T) 2 5X011 
VHASTROs TIDE s4Xe THSP(1I+1) p8X95HSP(I) eSXLIHWIND STRESS0/067X 0 6H(FFE 
2T) 0 OXO4H(FT) 0 LOX oe 4b ( FT) eX» SHPARAMETER?//) 
bo 390 LL=1leLMM 
READ (39) WWX(LL41) oWWX(LL) pWWY(LL41) oWWY(LL) o SAC(N) oSP(LL41) o SP(LL 
1) oWK(LL) P 
WRITE (60380) wwxX(LL4+1) ewwX (LL) owWwY(LL41) oWWY (LL) »SAC(N) oSP(LI +4) 0S 
IP(LL) pWK (LL) : 
3E0 FORMAT (1OX0F Tele 8XoF Tole TXoF Tele BX oF Toso TXoF4ol ol OXoFUo201l0KeF4.2 
1)55X0F10e3) 
390 CONTINUE 
460 CONTINUE 
405 CONTINUE 
REWIND 39 
WRITE (60410) 
IF (IOMIT .EG@e 1) GO TO 422 
410 FORMAT (iH1) 
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Cc 


Cc 


C 
Cc 


WRITE (69420) 
420 FORMAT (54XeLLHOUTPUT DATAe/*S4xolineerat2= e--=9///) 
WRITE (60260) POrePNeReWXe VE * SEs BFF ewKCOR? Cis C20C3 
422 CONTINUE 
DO 425 L=i»LMM 
INCREMENTAL X=VALUES AND AVGe DEPTHS ARE COMPUTED. 
DELX(L) = DIST(L) = DIST(L+1) 
DAVG(L) =(p(L) + D(Lt+1))/2.0 
425 CONTINUE 
DO 430 L=1leLM 
PHI(L) = (PHI(L) * PI)/180.U 
THE SINE OF THE LATITUDE IS EVALUATED. 
SINPHI(L) = SIN(PHI(L)) 
430 CONTINUE 
TII = 00 
DO 500 N=ieNm 
IF (IOMIT .E@e 1) GO TO 435 
WRITE (60410) 
435 CONTINUE 
READ (39) DELT(N) 
SUMSX = 020 
SUMSY = 0:0 
IF (IOMIT -EQ@e 1} GO TO 445 
Ties! Wert 
TII = T1 + DELT(N) 
WRITE (69360) TIeTII 
WRITE (60440) 


440 FORMAT (1Xe5HDIST.e 3X *SHDEPTH9 5x0 6Hi) AVGee 3X9 SHPRFSeRISE 0 2X0 1OHAST 
LROe TIDE? 3X0 VHINITIAL » SX OHYMFLUX 0 5X 9 THONSHORE 0 4X0 1OHALONGSHORE 0 2X 
210GHTOTAL WIND? 3X eSHTOTALe/o1X0 44: (NM) 0 3X9 BH(FT oMLW) 0 2X9 BH(FT eMLW) 0 & 
SXo4tH(FT) op 7XeGH(FT) poXe LOHLEVEL (ET) o LXeLLH(FT*FT/SEC) eo AXe 1 NHSFTUP 
YET) »2X010HSETUP (FT) e2xXe10HSETUP (FT) *2Xe10HSURGE (FT) 0//) 


445 CONTINUE 
DO 490 L=isotLmMM 
STORM SURGE COMPUTATIONS BEGINe 


READ (39) wWwx(L4+1) oWwwX(L) owWY (L+1) owWwY(L) »SA(N) oSP(L4+1) 9SP(L) oWK OL 


1) 

AN = WK(L) * (WWxX(L) + Wwx(L+1)) * WwKCOR 

EN = WK(L) * GWWY(L) + WWY(L41)) * 065 * wKCOR 
SPN = (SP(L) + SP(L4+1))/2-0 
IF (N eNEe 1) GO TO 450 ; 

PROBLEM IS IWNIALIZED FOR FIRST TIME LEVEL. 


SXP(L) = 0.60 
SYP(L) = 0.0 
SPP(L) = SPN 
VP(L) = 060 
BP(L) = BN 
SAP SA(N) 
450 OTS DAVG(L) + SE + SXP(L) + SYP(L) 


DTN = DTS + SA(N) + SPN 


bA = C3/(0TH) : 
DEN = DA *DA *BFF * ABS(VP(L)) *DELT (TH) 


VN = ((bN + BP(L))* DELT(N) * 0.5 + VP(L))/(1.0 + DEN) 


BY = ABS(BN) 
VNTEST = SQRT(BY/( BFF*DA*DA) ) 
VNN = ABS(VN) 


A CHECK IS MADE TO INSURE THAT THE ALONGSHORE NOSE NOT EXCEEN THE 


LIMITING VALUE 
IF (VNTEST eGTe VNN) GO TO 460 
IF (VN) 4560457,457 


50 


DTS + ((SA(N) + SAP)/2.U) + ((SPP(L) + SPN)/2.0) 


c 
Cc 


456 VN = > VNTEST 
6O TO 46) 
457 VN = VNTEST 
460 DELSYN = (C2 * NELX(L) -* (SINPHT(L) + SINPHI(L4+1)) * VN) /DTN 
DELSXN = (Cl * MELX(L) * AN) JOT 
SUMSX = SUMSX + GELSXN 
SUMSY = SUMSY + HELSYN 
VS = 774400 * VIN 
ST = SA(N) + SE + SUMS + SUMSY + SPN 
Sw = SUMSX + SUMNSY 
iF (IOMIT .EQ@e 1) GO TO 485 
TF. [OMIT NGT EQUAL Le GUTPUT RESULTS ARE PRINTED GUT IN NDETATL FOR 
EACH TIME LEVEL. 
WRITE (60470) DIST(L) eL(L) 
470 FORMAT (1XeF5SeleSXeF 501) 
wRITE (60480) DAVG(L) #SPN»eSA(N) 6 SE eVS0SUMSX oe SUMSY 6 SW ST 
YEO FURMAT (21K9FDeLeDKoF4 elo BX9FU4 ola TXoFu elo BXoF Tele SX0F6a50 TX 0FH0 306 


1X0F6030 5X0 50a) 
4&5 SXP(L) = SUMSX 
SYP(L) = SUMSY 
SPP(L) =, SPN 
VP(L) = VN 
BP(L) = AN 


400 CONTINUE s 
IF (IOMIT .EGe 1) GO TO 495 
WRITE (60470) DIST (LM) 0 D(LM) 
mEW VALUES ARE SET EQUAL TO PREVIOUS VALUES. 
495 Sx(N) = SUMSX 
SY(N) = SUMSY 
SwWw(Nn) = Sw 
SPRES(N) = SPN 
STOT(N) = ST 
500 CONTINUE 
Kena 
565 CONTINUE 
WRITE (6°410) 
WRITE (60510) 
THE RESULTS ARE SUMMARIZED FOR ALL TIMF LEVELS. 
510 FORMAT (38Xe41HSUMMAKY OF OPEN COAST STORM SURGE PROBLEM? ///) 
WRITE (601) (A(T) 0i=l°20) 
WRITE (69260) POoPNoReWXo VF OSE eiiFFowKCORPClieCeeC3 
WRITE (60520) 

520 FORMAT (9///02Xe4HTIME® 6X 0 SHSETUP 0 SX eo SHSETUP 0 4X0 9HTOTe WIND ® 4X0 GHA 
15Te TIDEs 3X1 SHINITIAL WATER ® 4X» SHPRESSURE o 3XeLLHTOTAIL WATFEFRe/ 01x09 
27H( HOURS) 9 3X0 THX =COMP © 0 3X0 THY=COMP 6 0 4X 9 QHSETUP » 7X 0 SHLF VEL 9 9X 9 SHLEV 
SEL 0 LOK 0 SHSETUP® TX eSHLEVEL 9/0 1L0Xe9H(F Te MLW) oLKoOH(F Te MLW) 0 2X e 9H (F 
4T. MLW) 0 3X09K(FT. MLW) o4xe9H(FT> MLw)e 7X09H(FT > MLw) e 3Xe9H( FT. wlwW 
5)0//) 

AHOUR = 0-0 

DO 540 N=lenM 

AHOUR = AHOUR + HELT(N) 

WRITE (69530) AHGUIKeSX(N) 6 SY Civ) oSWWON) eSACN) 0SE®SPRES(M) 0 STOT(N) 

5360 FORNAT (2X 0F60e2r POX0FOee 0X PF eH OK oF Delo TX OF Sole GX eF Solel NXeF Do 20 
L4X0F5e2) ; 

540 CONTINUE 

IF (K eEQe 2) GO TO 5duU 
K = Kt1l 
GU TO 505 
550 CONTINUE 
END 
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QUAST=TWwO=i)IMFINSTONAL, OPFN COAST STORM SURGE PROGRAM 


CHFSAPEAKE BAY ENTRANCE 


INPUT NATO 
“CENTRAL PRESSUrE = 27.57 IN. HGo PFRIPHFRAL PRFSSURF = 29-92 IN. HG. 
RADIuS TU MAXIMUM WIND = 35.0 N.M. MAXIMUM GRGOTENT WIND = 102060 MI/HR. 
TRANSLATION SPEED = 22.0 KNOTS THITIAL RISF = 2.50 FFFT 
BoTTOM FRICTION FaCTOR = .@025 WINO STRESS CORRECTION FACTOR = 1.10 
C1 = 203.0 c2 = 106.0 C3 = 52h0.n 


SEARED PROFILE OVER THE CONTIiWFNTAL SHFLF 


DISTANCE FROM OFPTH 

SHORE (AeMe) (FFFT) 
62.00 600.0 
60.00 350.0 
55.00 240.0 
50.00 178.0 
45.00 154.0 
40.00 130.0 
35.00 115.0 
30.00 102.0 
25200 4.0 
20.00 71.0 
15.0u 64.0 
10.00 56.0 
5.00 44.0 
3.50 40.0 
2.00 36.0 
1.00 33.0 
200 30.0 


HEM RE RAR eg RK eg RHR EK KD KK HEE ERK ERO ET ERE KAR KS RE KEE Ke KERR ET RETESET KER EEE EK EE 


TIKE (1) = 200 TIVE (T+ %ceLT) = 250 
Wwx(T+1) wik(I) wwY (141) wey (T) ASTRO. TIOF SP(1I+1) sett) WIND STRFSS 
(FEFT) (FT) (FT) PARAMFTFR 
2201 72608 36.7 3902 2.50 228 02h 10000110 
w21e1 “2561 30.8 3607 2250 228 22h 00000110 
170 2101 2504 3008 2250 227 02h 20000110 
14.2 717.5 20.0 2564 2250 227 227 20000110 
rile2 “14.2 lo.S 2006 2.50 226 227 00000110 
“Seo 711.2 12.4 16.3 2250 226 226 90000110 
003 “3806 Gel 12.4 2-50 226 228 200000110 
“God “605 ood 9el 2.50 225 22h 290000110 
208 ~4e4 4.1 603 2.50 225 225 0000110 
“leo ~2.8 2.3 4el 2.50 225 225 290000110 
--7 -1le6 1.0 203 2050 224 025 00000110 
-02 -.7 23 1.0 2.50 224 224 290000110 
mel 702 ol °3 2250 024 224 «00000110 
0 mel 20 el 2.50 224 224 00000110 
700 -.0 oY) of 2050 224 24 210000110 
°0 -.0 a) 20 2650 224 224 290000110 


RG ER ERK ERE GK RR RE K EERE EE REESE KERR AT ERE HERE EE ee ER EM EKER TREE EEK A EERE RS ERE EE EE 


TIME (T) = 050 TIME (T+HDELT) = 1.00 

WWXCT+1) wwACL) awy(1+1) wy (1) ASTRO. TINF SP(1I+1) SP(T) WINN STRESS 
; (FEET) (FT) (FT) PARAOMFTFR 
3500 3720 51.4 54.4 029 29 200000110 
3003 ~“35eu G44 51.4 02h 229 -00000110 
“2509 “3003 37.9 444 228 oar 200000110 
721.9 "25.9 31.9 37.9 228 22h 210000110 
1802 721.9 20-5 31.9 227 oan 00000110 
w14es “18.2 21.5 26.5 227 227 00000110 
711.8 “14.68 17.1 21.5 226 227 200000110 
“Del 711.8 13.2 1701 226 226 00000110 
“607 -9e1 9.8 13.2 226 22h 200000110 
“48 76.7 6.9 9.8 225 226 200000110 
“Bel “4.8 4.5 6.9 225 025 290000110 
“let “301 2.6 4.5 225 225 -00000110 
m1e4 71.8 2.0 206 225 225 200000110 
71.2 104 1.7 2.0 225 225 290000130 
Se) 1.2 1.3 1.7 224 225 210000110 
oeo7 09 1.0 1.3 02h rat 200000110 
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TIME (T) = 30050 TIME (T+DELT) = 31.00 
DISTe DEPTH D AVGe PRESeRISE ASTROeTINE INITIAL YorFLUX ONSHORE ALONGSHORF TOTAL WIND TOTAL 
(NM) (FToMLW) (FT oMLW) (FT) (FT) LEVEL(FT)» (FT*ET/SFC) SETUP (FT) SETUP (FT) SFTUP (FT) SURGE (FT) 
62.0 600.0 
- 475.0 225 2.50 250 1072A 2000 2007 2007- 3026 
60.0 350.0 
29500 025 2.50 250 2901 2000 2015 2016 3.27 
55.0 240.0 
20920 226 2.50 250 “9701 2001 -.022 -.021 3.24 
50.0 178.0 
166.0 226 2.50 250 “B00? 2002 -.061 -.059 3.20 
45.0 154.0 
142.0 226 2.50 250 “7102 2002 2.101 -.099 3.16 
40.0 130.0 a 
122.5 227 2.50 250 "6307 2003 -.142 -.139 3.13 
35.0 115.0 
108-5 227 2.50 250 “5b 2005 -.185 -.181 3.09 
30.0 102.0 
93.0 227 2.50 50 51.9 2006 =.229 -.223 3.05 
25.0 84.0 
7705 028 2.50 050 "4500 2008 -.275 2267 3.01 
20.0 71.0 
67-5 228 2.50 250 4007 e011 -.322 -.311 2.97 
15.0 64.0 
60.0 229 2.50 250 3707 2014 -.371 -.357 2.93 
10.0 56.0 
50.0 029 2.50 250 =320R 2018 ~.422 -.404 2.89 
5.0 44.0 
42.0 229 2.50 250 =?A0o5S 2020 -.438 e418 2-88 
3.5 40.0 
3600 030 2.50 250 2603 2022 2454 =.432 2.86 
2.0 30.0 
3465 730 2.50 250 2305 2023 464 e441 2.86 
1.0 33.0 
31.5 230 2.50 250 72005 2024 474 -.450 2.85 
A) 30.0 
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SUMMARY OF OPFN COAST STORM SURGE PROBLEM 


CHESAPEAKE BAY ENTRANCE 


CENTRAL PRESSURE = 27-57 INe HGe PFRIPHERaAL PRFSSURF = 29.92 INe HG. 
RabIUS 10 MaXIMUM KIND = 35.0 NoMe MAXIMUM GRANTENT WIND = 10260 MI/HR. 
TRANSLATION SFEEU = 22.0 KNOTS TNITIAI RISF = .50 FFET 
BOTTOM FRICTLON FACTOR = .0025 WIND STRFSS CORRECTION FACTOR = 1.10 
cl = 203.0 c2 = 106-6 c3 = 5280.0 
Tine SETUP SETUP TOle WIND ast. TIDE INITTAL WATFR PRESSURE TOTAL WATFR 
(HOURS) X=-COMP.e Y-COMPe SETUP. LEVFL LEVEL SETUP LEVFL 
(FT+s MLW) (FT. MLW) (FT. Miw) (FT. Miw) (FT. MLW) (FT. MiLW) (FT. MLW) 
250 -.u0 20U 250 250 224 3.24 
1eU0 -.00 200 2250 050 024 Bory 
1-50 -.06 200 2250 225 3.25 
2000 -.01 290 2.50 226 3025 
2°90 -.01 201 2050 227 32?6 
3eud ~eul 91 2050 228 o77 
3-50 -.U2 el 2250 228 3028 
4eU0 wee 2 2050 029 3.29 
4eS50 ~.03 22 2050 230 3.30 
5°00 ~.04 203 2050 032 3-31 
5250 -.05 204 2050 233 3.3? 
6°00 -.06 20S 2.50 2°34 3033 
650 a's) 007 2050 035 335 
7°00 10 209 2.50 037 3036 
7°50 7.13 el2 2050 039 3038 
&e00 -.16 e15 2050 240 3e40 
Be5U -.20 219 2050 242 342 
9eu0 7.224 224 2050 245 3045 
9e50 -.29 e351 250 247 3249 
10-00 -.55 235 205 250 3653 
10°50 ~o41 4B 2050 053 3.59 
11¢°0U -.49 259 2050 256 3266 
lledu ~.5b o72 2050 260 3.746 
12¢uU 2.04 286 2050 264 308A 
12°50 72 1.206 2050 269 4.06 
13°00 ~.74 1ed2 2050 075 432 
13°50 .b7 1.01 2250 2A? 4e76 
14euu 7-242 1.96 2050 290 Se44 
14°50 ell 2037 2250 099 6048 
15°vu0 1.0u 2.91 2.50 1.10 791 
15¢50 2.26 3.21 2050 1.94 971 
16°UU 5.75 3234 2050 1.39 11047 
1605u dell 3.uU5 250 1.55 12.074 
17-0U 6.09 2062 2.50 1.70 13041 
17¢5u Oe4h 1.99 2650 1.78 13.25 
18°06 6.19 Led4 2250 1.76 12.28 
18650 5.25 ode 2650 1.61 1068 
19¢0U 4.15 220 250 1646 BeAT 
19650 3.19 -.20 2050 1.31 7.25 
20¢UU 2.42 =. 72 2050 1017 50K 
e0e5u 1.45 -1.07 2.50 1.94 4eAS 
21°00 1.43 -1.30 2050 094 4oN6 
21°50 1.11 “1.44 2050 085 3052 
22°00 268 1.51 2.50 278 3.15 
22°50 e71 “1.52 2050 e71 2.90 
2300 055 1.48 2.50 °AG 2.76 
23050 047 1.42 2050 262 2067 
24-0u 259 “1.35 2-50 258 2062 
24250 232 -1.27 2250 054 2059 
25200 220 1.21 2050 251 2.057 
25°50 222 “1.14 2650 2048 2056 
26°00 old 71.06 2.50 046 205A 
26°50 els -.99 205U 043 2059 
27-00 e1l2 -.92 2050 041 200A? 
27°50 elu -.55 2050 05D 2064 
2aeUU U8 -.79 2650 238 2067 
2850 207 -e7e 2050 236 2.70 
29euU0 °uS -.66 2650 235 2074 
29050 204 -.61 2.50 233 2077 
30°u0 004 056 250 03? 2eA0 
30°5U 0U3 ~e51 250 e351 2083 
31°00 202 ~o47 205U 30 20A5 
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